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ABSTRACT 

Gravitational microlensing events with high peak magnifications provide a 
much enhanced sensitivity to the detection of planets around the lens star. How- 
ever, estimates of peak magnification during the early stages of an event by means 
of minimization frequently involve an overprediction, making observing cam- 



CN ■ paigns with strategies that rely on these predictions inefficient. 

^ ! I show that a rudimentary Bayesian formulation, incorporating the known sta- 

O , tistical characteristics of a detection system, produces much more accurate pre- 

Q ! dictions of peak magnification than minimisation. Implementation of this sys- 

^ ! tem will allow efficient follow-up observing programs that focus solely on events 

^! that contribute to planetary abundance statistics. 

O 

c/3 . Subject headings: gravitational lensing — methods: data analysis 

> 

ij ■ 1. Introduction 

P.. 

Following the suggestion of Paczyhski (1986), several collaborations, notably MACHO 
and EROS, began to search for gravitational microlensing towards the Magellanic Clouds as 
an indicator of compact objects in the halo of the Milky Way (Alcock et al. 1993; Aubourg et 
al. 1993). At about the same time, the OGLE collaboration began a survey in the direction 
of the Galactic bulge (Udalski et al. 1992, 1993). It was soon found that a much higher event 
rate occurred in fields towards the Galactic bulge relative to the rate towards the Magellanic 
Clouds (Udalski et al. 1994a; Alcock et al. 1995, 1997a). Since 1990, approximately 1000 
such events have been detected (Alcock et al. 2000; Udalski et al. 2000). 

Several groups including PLANET (Probing Lensing Anomalies NETwork, Albrow et al. 
1998, 2001; Dominik et al. 2002; Gaudi et al. 2002), MPS (Microlensing Planet Search, Rhie 
et al. 1999) and /xFUN (Microlensing Follow-Up Network, Yoo et al. 2004) monitor events 



2 



much more intensively than the survey groups in order to identify anomalous behavior that 
can signal the presence of a planet associated with the lens star. High-magnification events 
in particular (those with Aq > 10) attract the attention of follow-up groups since it is these 
that arc most likely to give detectable planetary signals (Griest & Safizadeh 1998; Gaudi, 
Nabcr & Sackctt 1998). In addition, for high magnification events the angular size of the 
source star may be non-negligible in comparison to the lens-source angular separation. In 
these cases the hghtcurves of the events can provide the possibility to determine the lens- 
source relative proper motion (Gould 1994; Alcock et al. 1997b) and atmospheric properties 
of the source (Heyrovsky 2003). 

In the first years of operation, when microlensing alerts came primarily from the MA- 
CHO collaboration, detected event rates were low enough that PLANET could monitor 
almost all potentially interesting events with ease. For the last two years (the 2002 and 2003 
Bulge seasons) , this has not been the case, due to the much improved alert rate since the ad- 
vent of the OGLE III early warning system (EWS), http://www.astrouw.edu.pl/ oglc/ogle3/ews/ews.html 
(Udalski et al. 1994b; Udalski 2003). In excess of 400 events were alerted by the EWS in each 
of these years. In addition, approximately 75 events were alerted in 2003 by the MOA col- 
laboration (Bond et al. 2002) although some of these were duplicates of EWS events. We are 
now in an era in which a careful selection of events is necessary to optimize planet detection 
and exclusion productivity. For this reason, follow-up groups require accurate predictions of 
eventual maximum amplications in the early days following a detection. For the remainder 
of this paper I will focus exclusively on events detected by the OGLE III EWS. 



2. Fitting microlensing lightcurves 



Most microlensing events are well fitted 
for the magnification A{t) at time t, 



by a point-source point-mass-lens (PSPL) model 



A{u) 



(1) 



u^/u'^ + 4' 




(2) 



r(t) 



t-to 



(3) 



The impact parameter u is the angular separation between the source and lens measured in 
units of the angular Einstein radius. 




(4) 
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where M is the mass of the lens, Ds is the observer-source distance, Dl the observer-lens 
distance, Dls the source-lens distance and Uq is the impact parameter at to, the time of 
maximum magnification. The Einstein radius crossing time, 

tE = —, (5) 

A'rel 

where n^ei is the relative proper motion between the lens and source. 

The lightcurvc of a PSPL event can thus be characterised by 3 + 2n parameters, 
{to,Uo,tE) plus for each n telescope -|- filter combinations, the unmagnificd (baseline) mag- 
nitude of the source star mbase and the blending parameter /bi, where 1 — fu is the fraction 
of blended (non-lensed) light. The maximum magnification, Aq — A{uo) frequently replaces 
uo as a parameter. 

The conventional method for predicting peak magnifications is to use minimization 
techniques to fit PSPL models to data from OGLE III (possibly supplemented by a follow- 
up group's own data) as they become available. Such fits are continuously updated and 
the subsequent predictions revised as data accumulate. Experience has shown that early 
predictions of eventual maximum magnification using these methods systematically yield 
overpredictions, strongly limiting the usefulness of such estimates. In particular, very large 
maximum magnifications (with large uncertainties) are often predicted for events that turn 
out to be of rather low amplitude. Valuable observing time is often wasted monitoring such 
events in order to confirm their nature. 

The reason for this overprediction of values for Aq is that in using minimization 
for a predictive purpose, one implicitly assumes that all parameter values are equally likely. 
However, for microlensing events this is far from being the case. From a purely geometrical 
perspective, high magnification events are exceedingly rare. In practice, being of high mag- 
nification, they have a higher probability of detection by a survey group. It is the individual 
detection efficiency of a survey convolved with the intrinsic event rates (both of these as a 
function of event parameters) that determines the magnification probability, given that an 
event has been detected. 

These ideas can be given a quantitative basis in a Bayesian formulation of the problem. 
The merits of the Bayesian approach to statistical analysis have been discussed at length 
elsewhere and will not be reargued (Loredo 1990; Sivia 1996). Here we simply note that a 
Bayesian formulation with appropriate priors should produce an unbiased estimate of the 
eventual microlensing event parameters during the rising part of a lightcurve. 

Prom Bayes' theorem, the probability density for a microlensing event to have a certain 
set of parameter values 6, given the fact, O, that it has been detected by the OGLE HI EWS 
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and that data D have been acquired, 

p{9\D,0)^piD\e,0)pie\0), (6) 

where ^ = {Aq, t^, to, ml^^, /^j, j = l..n). The second term on the right hand side of equation 
(6), p{0_\O) (known as the prior), is the underlying probabihty density for ^ given a detection, 
i.e. p{9\0)d9 is the probabihty that 9 is in the range [9,9 + d9]. It is this function that 
incorporates both the natural event occurrence probability and the particular parameter 
sensitivities of the detection system. The first term, p{D\9,0) — C{9\D), the hkelihood 
function for 9_ given D. Data from the OGLE III survey consist of /-band magnitudes and 
their uncertainties {mi,ai) at time ti (i.e. n — 1). I assume each rrii to be drawn randomly 
from a normal distribution N{mifi, (7i) where the true value of the magnitude at time ti is 
rriifi. This implies that 

C{9\D) oc e-^'/2 (7) 

where 



I mi-m(t„9) \ 
i \ / 

and m(ti, ff) is the magnitude evaluated from the model parameterised by 9_. 

Analagous to a minimization, the value of 9_ that maximises p{9\D,0) (i.e. the 
posterior mode) is taken as the best estimator of ^q, the true value of 9_. In the absence of a 
prior, this solution reduces to the minimum solution. We stress here that when sufficient 
data are available to constrain a fit to a certain event, e.g. when the event is over, the 
Bayesian and minimization techniques give the same parameter values and the choice of 
prior is largely irrelevant. In other words, the solution is not driven by prior probabilities 
when sufficent empirical information is available (see Sivia 1996, Chapter 2 for a discussion 
of this point). 

If the parameters 9_ are statistically independent quantities^ p{9\0) factorizes as 

p{9\0) = p{Ao) p{t^) p{to) p{m^) p(/bi), (9) 

where for brevity I have omitted " |0" in the probabihty densities on the right hand side of 
the equation. It is often more convenient to work in decadic logarithmic units for several of 
these quantities, in which case 

..|^^ p(lg Aq) p(lg(tE/r)) p(lg(Ato/r)) p(mb)p(/bi) 
"^-^^^ = (lnlO)3Ao (W^l (AtoAl ' ^ ^ 
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^This is not necessarily the case given that a detection system may preferentially select events with 
parameter corelations, however inspection of EWS-detected events has not revealed any such corelations as 
yet. 
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with t* being an arbitrary unit of time. Here I define Ato to be the time to peak magnification 
from an initial "alert date" . For the remainder of this paper I adopt t* — 1 d. 

It is worth noting that even if each parameter in 9 is independent in p{9\0), they are not 
independent in the likehhood function p{D\9,0) and hence not independent in p{9\D,0). 
Thus when fitting a model to a lightcurve, particularly when only early data are available, 
the fitted maximum magnification is affected not only by the prior on Aq but also by the 
priors on the other other parameters. 



2.1. Inclusion of a blending parameter 

The criterion used by the OGLE III EWS is that a blending parameter is only used 
when it is more than 3-a less than unity and is larger than its formal uncertainty. In this 
paper I use the odds ratio test, a natural way to decide between two different models. I 
define the odds ratio 

pmD,o) _ p{D\io)pmo) 

p{9\D,0) p{D\9,0)pi9\Oy ^ ' 

where 9_ = {Aq, Ie, to, mlg^^,j = l..n) indicates the set of model parameters without blending 
{n — 1 when considering only EWS data). Having no a priori indication about whether to 
include blending I choose p{9\0) — p{9\0). If we assume a unform prior probability density 
for /bl in the range < /bi < 1 and zero outside this range, and assuming a Gaussian 
probability density function for /bi about /bi,o, it can be shown (see for instance Sivia 1996 
Ch 4) that equation (11) reduces to 

Pm,0) ^ P{D\h) 1 .-.21 
p{9\D,0) p{D\9^) V2^af,, ^ ^ 

for cases in which /bi,o is more than several af^^ away from the cutoffs imposed by the prior. 
Otherwise, for /bi,o close to 1, equation (11) becomes 



p(g|D,0) p(D%) 1 



P(e\D,o) + 

while for /bi,o close to 

p{9\D,0) _ p{D\~9,) 1 

pi9_\D,0) Mi^l^o)^.,,,(l + erf(^)) 



(13) 



(14) 



The odds ratio is then made up of two terms. The first of these represents a relative "goodness 
of fit" between the two models while the second is the "Occam penalty" for introducing a 
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new parameter. Only when /bi,o < 1 and the odds ratio is less than unity is a blending 
parameter used in this analysis. 



I have used the set of microlensing events detected by the EWS in 2002 to determine 
the parameter priors. Since our interest is in the set of PSPL events, I have removed 41 
events that showed deviations from PSPL behavior from this analysis. Excluded events were 

numbered 18, 23, 40, 51, 68, 69, 77, 80, 81, 99, 113, 119, 126, 127, 128, 129, 131, 135, 143, 
149, 159, 175, 194, 202, 203, 205, 215, 228, 229, 232, 238, 254, 255, 256, 266, 273, 307, 315, 
339, 348, 360, out of the complete set of 389 alerts. 

Parameter values have been obtained for these events using a simplex downhill method 
to minimise (EWS estimates of Aq, tg, mt, fu can also be obtained from the EWS web 
page). For Ato we require an objective definition of an "alert date" that can be applied to 
all events. 1 have arbitrarily chosen a working definition of an alert date as being the date at 
which three successive data points have been more than 1-a brighter than m,],. the baseline 
magnitude. In practice, rrib can be determined separately from and in advance of the other 
parameters. 

The distributions of IgAo, Ig^E and IgAto are shown in Figure 1. For the purposes 
of obtaining Bayesian prior probability densities for these quantities, the distribution func- 
tions are adequately represented by the following empirically-chosen functions, also shown 
in Figure 1: 



It is also instructive to examine the distribution of Uo, shown in Figure 2(a). In the 
absence of any selection effects, this distribution should be uniform. In fact, there is an 
enhanced sensitivity to detection of high magnification (low Uq) events and a rapid decrease 
in sensitivity for uq > 0.85. Figure 2(a) also shows the shape of the adopted prior on Ig 
(eq. 15) when transformed to Uq. 

Figure 2(b) shows the same data but excluding those events where Aq has a formal 
uncertainty greater than 50%. This illustrates that many high amplification events have 
maxima that are poorly constrained from OGLE data alone. 



3. Statistical properties of the 2002 OGLE events 



PiUtE/t*)) 

p{lg{Ato/t*)) 



0.660 exp [-1.289 Ig Aq] 

0.476 exp [-{lg{tE/t*) - 1.333)70.330] 

0.156 exp [-(\g(Ato/t*) - 1.432)70.458] . 



(15) 
(16) 
(17) 
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4. Application to 2003 OGLE alerts 

As a test of the Bayesian method, I have apphed a fitting procedure that maximises 
p{9\D, O) to a sample of the PSPL events alerted in real time by the OGLE III EWS in 2003. 
These consist of events OGLE-2003-BUL-138 to OGLE-2003-BUL-462 and excluding events 
numbered 145, 160, 168, 170, 176, 192, 200, 230, 236, 252, 260, 266, 267, 271, 282, 286, 293, 
303, 306, 311, 359, 380, 419 that do not appear to be due to PSPL microlensing and 188, 197, 
245, 263, 274, 297, 387, 399, 407, 412, 413, 417, 420, 422, 429, 430, 432, 433, 435, 437, 440, 
441, 442, 443, 444, 449, 450, 452, 453, 454, 455, 457, 459, 461, 462 that were still ongoing at 
the time of writing. Events OGLE-2003-BUL-137 and earlier were anounced by the EWS in 
a single email at the beginning of the 2003 Bulge season and thus not alerted in real time. 
OGLE-2003-BUL-238 (A. Gould 2004, private communication) and 262 (Yoo et al. 2004) 
are events in which the lens is known to have transited the source and OGLE-2003-BUL-208 
and 222 may also involve finite source effects. These events have not been excluded. For the 
remaining sample of 267 events, I have used only the OGLE III data taken before the EWS 
alert time, defined as the reception of the alert email by the author. For the zero point of 
A^o for each event, I have used the definition in § 2 except for cases in which this has not 
occured before the EWS alert time in which case the latter has been used as the zero point. 

As mentioned in § 1, different fitting codes can produce different estimates of maxi- 
mum magnifications, particularly for high-magnification events for which blending may be 
involved. In particular, there is a concern that a direct comparison of predictions with the 
EWS alert predictions may suffer from such differences. In order to compare the maximum 
magnifications predicted by the Bayesian method with those predicted using fitting, I 
have thus used very similar computer codes to make p{0_\D,O) and optimisations, elect- 
ing not to use the EWS-fitted parameters. To avoid the problem of slightly different blending 
parameters resulting in large differences in derived magnifications, for each event I compute 
the brightness increase. Am = mbase ~ ''^(^o); where m(to) is the magnitude at tg. The pre- 
dicted values of Am using only the pre-alert data for an event are compared with the values 
determined using all the data. When all the data are available, the parameters derived using 
and p{9\D, O) are almost always identical. Exceptions to this are in a few cases for which 
there are no data over the peak to constrain the fits. 



4.1. Comparison of Bayesian and predictions 

The predictive performances of the Bayesian and models at the time of EWS alert 
are illustrated in Figures 3-5. Figure 3 shows the distributions of predicted peak magni- 
fications for both models and compares these with the eventual values. Figure 4 shows the 
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same data as a function of uq. For the models, there is clearly a population of low uq 
events with predicted brightenings of more than 10 magnitudes that do not eventuate. Such 
overpredictions are not present in the Bayesian fitted models. On the other hand, there is 
a tendency for the Bayesian models to underpredict the peak, and at alert time to fail to 
predict the small population of high magnification events in Figure 3(a). In Figure 5 1 com- 
pare the distribution of the differences in predicted vs actual brightenings for both models. 
Again, the tendency for the fits to overestimate the peak is obvious. 

5. Case studies 

As pointed out in § 4.1, Bayscsian solutions to early lightcurve data often fail to indicate 
the nature of high magnification events. It would be of concern if high magnification events 
were not observed due to this tendency. To illustrate in more detail the behavior of Bayesian 
vs models, I consider here examples of low and high magnification events. These examples 
show several generic aspects of how Bayesian vs solutions evolve as data accumulate. 

OGLE-2003-BLG-171 was a low magnification event (Aq = 1.37). At this magnification, 
the source star barely passes within one Einstein radius of the lens and the event is unsuitable 
for detecting a planetary anomaly. This is typical of the type of event that a follow-up 
program should avoid observing. Figure 6 shows Bayesian fitted lightcurves at 5 day 

intervals as the event evolves from its alert date. The predicted maximum magnifications 
corresponding to each panel arc listed in Table 1. At alert, the event is predicted to be of 
low magnification but by JD 2452785 (panel d in Fig. 6), the solution suggests a high 
magnification, albeit with a large uncertainty. Since the lightcurve appears to be rising 
rapidly, follow-up programs may well begin observing the event in order to improve on 
the high uncertainty in the predicted peak magnification. As more data accumulate, the 
low-magnification nature of the event becomes apparent. Although the true nature of the 
event would be identified relatively quickly by a follow-up observing program, there is a 
not-insignificant overhead associated with adding the event to the program. In contrast to 
the x^ fit, the predicted peak magnification for the Bayesian solution changes steadily with 
time. At no time is a high-magnification event suggested and a follow-up strategy based on 
this method would ignore the event. 

OGLE-2003-BLG-208 (Fig. 7, Table 2) reached a moderately high magnification {Aq ^ 
45, Ao, unblended ~ 17). Thc projected source trajectory passed as close as 0.02 9e to the lens 
and thus had a high probability of intersecting a central caustic if it were present. Thc alert 
date for this event corresponds to panel (c) in Figure 7 at which time the predictions of 
peak magnification are 2 x 10^ and 4.4 for the x^ Bayesian solutions respectively. The 
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Bayesian prediction of Aq — 4.4 is sufficiently fiigh to warrant the attention of a follow-up 

observing program such as PLANET. As data acumulate, the predicted peak magnifica- 
tion rises until reaching 5 x 10^ in panel (g). The true peak magnification {Aq ~ 48 starts 
to become apparent from panel (h) as the event peaks. In contrast, the Bayesian predicted 
peak magnification rises steadily until the true peak magnification is identified from around 
the time of panels (f) - (g). 

The behavior illustrated by these two examples is typical. For low magnification events, 
the Bayesian model never indicates them as being worthy of observational follow-up. For 
high magnification events that should be observed, the peak magnification is initially un- 
derestimated but adjusts to an appropriate prediction as soon as the data indicate. In all 
cases examined, this occurs relatively early in the event when the magnification, A <3. For 
both high and low magnification events, the Bayesian predicted peak magnification changes 
smoothly while the prediction is prone to large changes as new data points are included. 
The Bayesian solutions usually converge to the correct amplification earlier than the 
solutions. 



6. Summary 

High magnification events provide the best opportunity for detecting signals of planets 
around lens stars and for obtaining upper limits on their abundances. Intensive photometric 
monitoring programs arc hampered currently by difficTiltics in identifying high magnification 
events well before peak. Systems that use minimization to fit PSPL models to early data 
are prone to exagerated predictions of peak magnification. Such prections induce observers 
to spend their time monitoring events that ultimately have little statistical power. 

I have shown here that a predictive system based on a Bayesian formalism that takes 
account of the characteristics of a detection system is immune to such behavior. Although 
such a Bayesian system tends to initially underpredict the peak for high magnification events, 
accurate prediction occurs as soon as sufficient data accumulate to justify the assertion. In 
all cases examined, this occurs well ahead of peak in their associated lightcurves and early 
enough for the events to be targeted for observation. Implementation of such a system based 
on the OGLE Early Warning System should result in much improved observing productivity 
for the 2004 season. 

I am grateful to Martin Dominik for his comments on an earlier version of this paper. 
I think the referee, Andy Gould, for his suggested improvements to the manuscript. This 
work was supported by the Marsden Fund under contract UOC302. 
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Fig. 1. — Distributions of (a) peak magnifications (b) Einstein timescales and (c) times from 
alert to peak magnification from my own fits to the 2002 OGLE event data. The sohd fines 
are the adopted Bayesian priors based on these distributions. 
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Fig. 2. — Fractional distribution of uq for the 2002 OGLE PSPL events. Panel (a) includes 
the entire data set while panel (b) excludes events for which the uncertainty in Aq is greater 
than 50%. The solid line in panel (a) represents the Bayesian prior on Ig^o transformed to 
uq and scaled to the first data point. 
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Fig. 3. — Fractional distributions of (a) OGLE peak magnifications with (b) predictions 
and (c) Bayesian predictions at the time of alert. In panel (a) the solid line represents the 
fits and the (mostly overplotted) dotted fine represents the Bayesian fits. 
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Fig. 4. — Fractional distributions of Uq for the 2003 OGLE PSPL events. Panel (a) shows 
the values derived using all the data while panels (b) and (c) are respectively the ^-iid 
Bayesian predictions at the time of alert. 



-16- 



0.3 
0.2 








1 















3 


:ion 





2 


o 






o 





1 












.0 



10 

m 





m (t 



1 

J predicted 



,0 



Fig. 5. — Fractional distributions of predicted peak magnitudes minus eventual peak mag- 
nitudes for (a) and (b) Bayesian models. 

Table 1. Peak magnifications for the and Bayesian fits to OGLE-2003-BLG-171 

corresponding to the panels in Figure 4. 



Bayesian 



Panel 


A 




/bl 




■^unblended 


A 




/bl 




^unblended 


a 


1.240 


0.670 


1.000 




1.240 


1.053 


0.014 


1.000 




1.053 


b 


1.050 


0.009 


1.000 




1.050 


1.051 


0.009 


1.000 




1.051 


c 


1.050 


0.009 


1.000 




1.050 


1.051 


0.009 


1.000 




1.051 


d 


34934.748 




1.000 




34934.748 


1.300 


0.314 


1.000 




1.300 


e 


1.449 


0.557 


1.000 




1.449 


1.273 


0.150 


1.000 




1.273 


f 


1.459 


0.330 


1.000 




1.459 


1.367 


0.169 


1.000 




1.367 


g 


1.539 


0.353 


1.000 




1.539 


1.446 


0.199 


1.000 




1.446 


h 


1.404 


0.082 


1.000 




1.404 


1.392 


0.072 


1.000 




1.392 


i 


2.524 


1.366 


0.233 


0.206 


1.355 


1.356 


0.019 


1.000 




1.356 


j 


1.374 


0.005 


1.000 




1.374 


1.374 


0.005 
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1.374 




Fig. 6.— Evolution of Bayesian (solid line) and (dashed line) fits to OGLE-2003-BLG- 
171. In panels (b), (c), (i) and (j) the dashed line is overprinted by the sohd fine. Axis 
ranges are the same for all panels. 




Fig. 7.— Evolution of Bayesian (solid line) and ^ (dashed line) fits to OGLE-2003-BLG- 
208. In panels (h) - (j) the dashed fine is overprinted by the sohd fine. Axis ranges are the 
same for all panels. 
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Table 2. Peak magnifications, blend fractions and deblended peak magnifications for tlie 
and Bayesian fits to OGLE-2003-BLG-208 corresponding to tlie panels in Figure 5. 



Bayesian 



Panel 


A 




/bl 




■^unblended 


A 




/bl 




■^unblended 


a 


63601.415 




1.000 




63601.415 


2.055 


0.657 


1.000 




2.055 


b 


138528.609 




0.209 


0.349 


28987.904 


2.150 


0.700 


1.000 




2.150 
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